Temperature of the inflaton and duration of inflation from WMAP data 
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If the initial state of the inflaton field is taken to have a thermal distribution instead of the 
conventional zero particle vacuum state then the curvature power spectrum gets modified by a tem- 
perature dependent factor such that the fluctuation spectrum of the microwave background radiation 
is enhanced at larger angles. We compare this modified cosmic microwave background spectrum 
with Wilkinson microwave anisotropy probe data to obtain an upper bound on the temperature 
of the inflaton at the time our current horizon crossed the horizon during inflation. We further 
conclude that there must be additional e-foldings of inflation beyond what is needed to solve the 
horizon problem. 



The inflationary paradigm [1] provides a ready solu- 
tion for the horizon and flatness problems while simul- 
taneously providing a source for the initial density per- 
turbations required as a seed for the large scale structure 
in the universe today. During inflation the energy den- 
sity of the universe is dominated by the potential energy 
of a scalar field leading to a brief period of accelerated 
expansion. 

The amplitude for quantum fluctuations in the inflaton 
field <p generated during inflation effectively freezes out 
for a particular comoving scale once that mode crosses 
the horizon. These fluctuations also translate into per- 
turbations in energy density and curvature in the uni- 
verse. In the inflationary paradigm these perturbations 
are the source of the observed anisotropy in the mi- 
crowave background and the seed for large scale struc- 
ture. Conventionally the fluctuation power spectrum for 
the inflaton is calculated using the zero particle vacuum. 
However if inflation was preceded by a radiation era then 
the inflaton was in thermal equilibrium at some point 
in the early universe, via at least gravitational couplings 
at the Planck scale. In chaotic [2] and natural inflation 
models with a symmetry breaking scale / ~ Mpi [3] , non- 
zero momentum modes can be in thermal equilibrium due 
to gravitational interactions with the existing radiation. 
The inflaton (or any other field which contributes to the 
density perturbation during inflation) can then have a 
thermal distribution corresponding to the temperature % 
at the beginning of inflation. Though the inflaton field 
may be decoupled from the rest of radiation before infla- 
tion it retains its thermal distribution in an adiabatically 
expanding universe. (This is different from the warm in- 
flation scenarios [4] where the inflaton decays to radiation 
during inflation.) In reference [5] it was argued that the 
power spectrum of gravitational waves generated during 
inflation has an extra temperature dependent factor due 
to the thermal gravitons which decouple at the Planck 
era and retain their thermal distribution. 

In this Letter we include the initial thermal nature of 
the inflaton <f> in the generation of the inflaton fluctua- 
tion power spectrum. The power spectrum of the pertur- 
bations of the inflaton, and hence of the curvature, has 



an extra temperature dependent factor coth[fc/(2<Zj7i)] 
where k is the wavenumber of the modes in comoving 
coordinates and a, is the scale factor at ij when infla- 
tion commences. The physical wavenumber is k/a and 
physical temperature T = T/a, where T is the comoving 
temperature during inflation. Using this thermal power 
spectrum we calculate the CMB anisotropy using CMB- 
FAST [6]. We compare the result with data from WMAP 
[7] and find that the comoving temperature of the infla- 
ton is constrained to be T < 1.0 x 10 _3 Mpc~ . Since 
T = 7o a Oj where Tq and ao are the temperature and 
scale factor when our current horizon scale crossed the 
horizon during inflation, this constraint can be rewrit- 
ten as Tq < 4.2H, where H is the Hubble parameter 
during inflation (all bounds stated in this paper are at 
ler or 68%C.L.). This result is valid independent of the 
scale of inflation. If inflation takes place in the GUT era 
(M GUT ~ 10 15 GeV) then % < 1.0 x 10 12 GeV and for 
the inflation at the electroweak era (Mew ~ 10 3 GeV) 
[8] we have the bound T < 1.0 x 10" 12 GcV. Since T 
is much less than the energy scale of inflation, Mi n f, 
the duration of inflation must be longer than what is 
needed to solve the horizon problem in models where the 
temperature at the beginning of inflation % ~ Mj„^/7, 
the energy scale of inflation. If AN is the number of 
e-foldings from the beginning of inflation to the time 
our current horizon scale crossed the de Sitter hori- 
zon, then the bound on Tq from WMAP data implies 
AJV = \n{Ti/%) > ln(0.03M in/ /ff). For GUT scale 
(M in f = 10 15 GeV) inflation AN GUT > 7 and for infla- 
tion at the the electroweak scale ANew > 32. 

We now calculate the curvature power spectrum due 
to inflaton fluctuations. Presuming a quasi de Sitter uni- 
verse during inflation, conformal time r (dr = dt/a) and 
the scale factor during inflation a(r) are related by a(r) = 

-lAHYfl - el where e = = = Mh ( ¥L) 

V is the inflaton potential, V stands for the derivative of 
the potential with respect to <p and H implies a derivative 
of H with respect to cosmic time. Henceforth a dot will 
always signify a derivative with respect to cosmic time 
whereas a prime may mean a derivative with respect to 
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the conformal time when it appears with quantities which 
are explicit functions of time or a derivative with respect 
to the inflaton field when it appears with the potential. 
Mpi is the Plank mass. We define additional slow roll 



parameters as rj = {^rj an( i $ = 4>/{H4>) = e — r/. 

We further define z = and u = a(r)<50(x, r), where 
<5</>(x, t) is the inflaton field perturbations on spatially flat 
hypersurfaces. Then the action for u is given as [9-11]: 



S = 



drd^x 



(«) 



/\2 



(Vii) 2 + — u 2 

z 



(1) 



The gauge invariant comoving curvature perturbation is 
defined as TZ = tb — H$$- where ib is the usual scalar 

metric perturbation. Consequently on spatially flat hy- 
persurfaces u = —zlZ. Expressing w(x, r) as a quantum 
field we can write: 



r) = J ^ (a k / fe (r) + a f _ k 



(2) 



where k is the comoving wavenumber. 

The gauge invariant comoving curvature perturbation 
can then be expressed as: 



z 



d 3 k 



(2tt) 3 / 2 
d 3 k 

^ k (r)e lkx 



fk(r)) e* k ' : 



The power spectrum of the comoving curvature pertur- 
bations can be defined by the relation 



2-7T 2 



k') 



(4) 



The usual quantization condition between the fields and 
their canonical momenta yields [ak, a£,] = S 3 (k — k') and 
the vacuum satisfies ak|0) = 0. If the inflaton field had 
zero occupation prior to inflation then (a^at) = and we 
would obtain a correlation function ~ |,A;(t)| 2 . However 
if the inflaton field was in thermal equilibrium at some 
earlier epoch [12] it will retain its thermal distribution 
even after decoupling from the other radiation fields and 
its occupation number will be given by: 



(4a k >) 
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<5 3 (k-k'). 



(5) 



where Ek is the energy corresponding to the k mode at 
the inflaton decoupling temperature Tf. For effectively 
free modes Ek ~ k/aj, and so Ek/Tf = k/(ajTf) = 
k/{a i T i ) = k/T. 

Using Eq. (3) and Eq. (5) it can be seen that 
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|/ fe (r)| 2 <5 3 (k-k') 



k 

2T 



5 3 (k-k>) (6) 



From the defining relation Eq. (4) for the curvature power 
spectrum and Eq. (6) we find that the power spectrum 
for the curvature perturbations can be expressed in terms 
of the mode functions /fe(r) as 



Pn(k) 



k 3 \h\ 2 
2tt 2 |z| 2 



coth 



2T 



(7) 



For constant e and S the mode functions /fe(r) obey the 
minimally coupled Klein-Gordon equation [11,13], 



k - ~ 

T 



fk = 0, 



where k = |k| and, for small e and S, v 
Eq. (8) has the general solution given by, 



(8) 



+ 2e + S. 



fk(r) 



Cl 



(k) HpX-kr) + C2 (fe) Hi 2 \-kr) 



(9) 



When the modes are well within the horizon they can 
be approximated by flat spacetime solutions /fe (t) = 
-^=e _lfcr , (k 3> aH) . Matching the general solution 
in Eq. (9) with the solution in the high frequency ("flat 
spacetime" ) limit gives the value of the constants of in- 



(3) tegration ci(fc) 



sf2L p i(y+\)\ 



and c 2 (k) = 0. Eq. (9) 



then implies that for — fcr » 1 or k <gC aH, 



(10) 



Substituting the solution in Eq. (10) for the the super- 
horizon modes (k <C aH) in the expression Eq. (7) for 
the curvature power spectrum we obtain 



Pn(k) 



H 4 



4tt 2 2 



/ k N 


n s -l 


' k ' 




| coth 




[aH, 


2T 



(11) 



with n s — 1 = 3 — 2^ = — 4e — 28. Note that since 
/fe(r) ~ z(t) for superhorizon modes, P-jz(k) is time in- 
dependent. 

We can now rewrite the power spectrum as, 



P n (k) - A(k ) 



/fc N 


Hs-l 


' k ' 




J coth 


\k 0/ 




2T 



(12) 



where ko is referred to as the pivot point and A(ko) is 
the normalisation constant. 



A(k ) 



H 4 



47T 2 ; 



; 2 \aH=k 



(13) 



where Hk a is the Hubble parameter evaluated when 
aH = ko during inflation. 
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During inflation a mode k leaves the horizon when 
k/a = H. Therefore for any k the initial ther- 
mal effects modify the scale-free power spectrum by 
coth(ofc-fffc/2T) = coth(iJfe/27fe), where ak and Tk are 
the scale factor and temperature when the mode k leaves 
the horizon. 

We now use the curvature power spectrum as given in 
Eq. (12) to obtain the CMB fluctuations. We make use 
of the publicly available code CMBFAST [6] and modify 
the power spectrum formula, according to our need, in 
subroutine powcrsflat in the program cmbflat.F. We set 
the parameters ilb — 0.046, fl m = 0.27, h = 71.0 as de- 
termined in [14] for a flat universe with no running of the 
scalar index n s . We set fc = 0.05 Mpc - . 

Fig. 1 shows two different possible CMB TT- 
correlation curves for two comoving inflaton tempera- 
tures, with n s = 0.99 and A(k ) = 2.01 x 10~ 9 . It shows 
clearly that as the comoving temperature of the infla- 
tons decreases we get a better fit of our TT-correlations 
with the WMAP data. The curves show that the change 
of temperature T of the inflatons essentially affects the 
lower multipoles or higher length scales of the TT- 
correlations while the higher multipoles remains unaf- 
fected. This is expected because coth(x) ~ 1 for x 3> 1. 

We now do a 3-parameter analysis by varying A(k ) 
and T for multiple values of n s from n s = 0.90 to 
n s = 1.1, and each time we compare the resultant TT- 
anisotropy with the WMAP data [7] and calculated the 
X 2 of the resultant computed distribution with respect to 
the WMAP data. From this three parameter fit we find 
that x 2 minimum occurs at the values A(ko) = 2.0 x 10~ 9 , 
n s = 0.99 and at the comoving temperature of the infla- 
ton field T = 0. In Fig. 2 we show the lcr allowed 
regions of T and A(ko) at different values of n s . The re- 
gions inside the curves are allowed by WMAP data, with 
the best fit value being A(k ) = 2.0 x 10~ 9 , n s = 0.99, 
and comoving temperature T = 0. The maximum al- 
lowed values for T increases when we increase n s from 
n s = 0.90 to n s = 1.06. For values of n s greater than 
1.06 maximum allowed value of T decreases monotoni- 
cally. So from this three parameter fits we see that the 
upper bound on the comoving temperature of the inflaton 
field is T < 0.001 Mpc" 1 . 

The power spectrum of tensor perturbations also has 
an identical coth(fc/2T) factor due to the thermal spec- 
trum of the gravitons in the initial state [5] . The tensor- 
to-scalar amplitude ratio r = 16e < 1.28 (95%C.L.) [15]. 
In our parameter fits we do not violate this upper bound 
on e which is independent of temperature. 

The temperature of the inflaton (Tq) when our present 
horizon was leaving the de Sitter horizon (H^ 1 ) can be 
calculated, T = %aQ < 4.2i? ? ~ 1 (where ao is the scale 
factor at the time of this horizon crossing). We then 
obtain 



To < 4.2 



&now \ D — 1 

a 



4.2H. 



(14) 



where a now is the current scale factor and Rh/dnow = 
H~ l /a . H = (87t/3) 1 / 2 M2 i/ /M p( , ignoring the vari- 
ation in H during inflation. The constraint Eq. (14) 
is valid whatever be the scale of inflation. If infla- 
tion takes place in the GUT era {M GUT - 10 15 GeV) 
then Tq < 1.0 x 10 12 GeV, and for inflation at the elec- 
troweak scale {Mew ~ 10 3 GeV) we have the bound 7o < 
1.0 x 10~ 12 GeV. In general, for M inf < 1.4 x 10 17 GeV, 
T <%. 
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FIG. 1. The dotted points with the error bars corresponds to 
the WMAP binned data for the CMB TT-anisotropies. The 
solid and dashed curves are generated from CMBFAST with 
the inflaton comoving temperature T = (2500Mpc) _1 and 
T = (25000Mpc) _1 respectively. The error bars on the WMAP 
data consists of measurement errors and errors attributed to cos- 
mic variance. 

From Fig. 2 we see that T has an upper bound at la- 
given by 

T < 1.0 x 10~ 3 Mpc" 1 = 4.2.R- 1 (68.3%C.L.) (15) 

where we have expressed the bound on the comoving 
temperature in terms of the present horizon scale Rh — 
4200 Mpc. 

The minimum number of c-foldings of inflation re- 
quired to solve the horizon problem is obtained by taking 
our current horizon scale to be the first to cross the hori- 
zon during inflation. If there is a thermal background of 
the inflaton and inflation commences at % ~ M in j/7 our 
results above imply that our current horizon must have 
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crossed the de Sitter horizon some time after the begin- 
ning of inflation. Thus the duration of inflation must 
be longer than what is just needed to solve the horizon 
problem by AiV = 
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FIG. 2. Values of T and A(ko) allowed at la by WMAP data 
are the regions inside the curves. The plots at different values of 
n s show that the maximum allowed value of T is 0.001 Mpc -1 . 



The upper bound on T from WMAP then implies 



AN > In 



0.03- 



M, 



inf 



H 



In 



0.01- 



M 



pi 



inf 



(16) 



Perturbations generated during these additional c- 
foldings are yet to enter our present horizon. For GUT 
scale inflation WMAP data implies ANqut > 7 and for 
inflation at the the electroweak scale it implies ANew > 
32. For Planck scale inflation, such as chaotic and natu- 
ral inflation, no additional e-foldings of inflation are re- 
quired. 

There exists a Lyth bound [16,17] on the minimum 
variation in the inflaton field during inflation, namely, 
5(f) > Mpi y/r/l SN, where 5 N is taken to be the num- 
ber of c-foldings of inflation (<~ 4.6) over which modes 
corresponding to rnultipolcSj I < 100 leave the horizon. 
This bound is further strengthened by a factor of 2.5 by 
including the additional 7 e-foldings (for GUT scale infla- 
tion) required by our analysis. The revision of the Lyth 
analysis by including the entire duration of inflation [18] 



should also include the additional c-foldings given in Eq. 
(16). 

In conclusion, we find that if inflation is preceded by 
a hot radiation era then there will be a stimulated emis- 
sion of inflaton perturbations during inflation into the 
initial thermal bath of the inflaton which will change 
the scale free spectrum. The change in the spectrum 
is large at large angles and there is no change in the 
CMB anisotropics at small angles. (This differs from the 
warm inflation scenario where the departure from scale 
invariance is proportional to the derivative of the infla- 
ton dissipative term, which is small to meet the slow roll 
condition [19].) Furthermore, from the WMAP data we 
find that the minimum number of c-foldings of inflation 
is larger than that required to solve the horizon problem, 
with the increase varying from 7-32 for the energy scale of 
inflation varying from the GUT to the electroweak scale. 
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